Refocusing schemes for holonomic quantum computation in presence of dissipation 
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The effects of dissipation on a holonomic quantum computation scheme are analyzed within 
the quantum-jump approach. We extend to the non-Abelian case the refocusing strategies formerly 
introduced for (Abelian) geometric computation. We show how double loop symmetrization schemes 
allow one to get rid of the unwanted influence of dissipation in the no-jump trajectory. 

PACS numbers: 03.67.Lx, 03.65. Vf, 03.65.Yz 



o 
o 

(N 



> : 

m . 
■ 

cn : 
o . 
^t- ■ 
o ■ 

^ : 

Ok 
-I— > ' 
G ■ 

a ; 

s : 



Due to the fragility of quantum coherence in the pres- 
ence of noise, decoherence is the main obstacle to the 
practical realization of quantum computation. Explor- 
ing potential ways to implement robust quantum com- 
putation therefore is a crucial and attractive challenge. 
To the aim of stabilizing quantum information a vari- 
ety of decoherence-reduction techniques have been devel- 
oped, such as quantum error correcting codes [H Q 
decoherence- free subspaces 0, and generalized decou- 
pling techniques @. 

The geometric/holonomic approach to quantum com- 
putation is believed to provide an intriguing geomet- 
rical way to protect quantum information processing 
0, H ES [H [H Geometric quantum effects, 

either the simple Abelian phase factor or its gen- 
eral non-Abelian counterpart |15| , arise when a quantum 
system is adiabatically driven to undergo an appropriate 
cyclic evolution. A purely geometric transformation can 
be obtained via a refocusing scheme @ which cancels out 
the dynamical phase or by properly designing the system 
levels 0, t° avoid the dynamical phase directly. 

The robustness of geometric quantum computation has 
been intensively addressed recently 0, 0, 0, 0, |2(], 
I2H ]. For a classical noise, the inherent stability of Berry 
phase has been clearly demonstrated by studying a spin 
one-half system subject to a randomly fluctuating mag- 
netic field [13 . 

In this paper we will consider geometric manipulations, 
Abelian and non Abelian, in presence of quantum dissipa- 
tion. We will first analyze, the quantum-jump formalism, 
the non-unitary effects induced by dissipation on the ge- 
ometrical gates in order to devise strategies to get rid of 
them. We will show that the geometric gate operation 
based on a double-loop refocusing scheme does possess a 
certain resilience against dissipation induced distortion. 
The main contribution of this paper will consist of a novel 
scheme which further extends the refocusing technique to 
enact non-Abelian holonomic quantum gates by means of 
a sort of simultaneous double-loop control process. 

Influence of dissipation on geometric quantum compu- 
tation. Let us consider the dissipative dynamics of a 
time-dependent periodic system governed by the follow- 
ing master equation 
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p = -i[H(t), P ] - - ]T(rtrv + P r{r k - 2T kP rl). (i) 

fe=l 



In the framework of quantum jump approach |22|. the 
operator T k is associated with jump in the trajectory; the 
time evolution of the system is calculated by averaging 
over all the possible trajectories. Stochastic jumps and 
diffusion of trajectories will destroy quantum coherence. 
This undesired effect becomes particularly severe for slow 
adiabatic evolutions since decoherence has a long time to 
bite in the process. 

As one could intuitively anticipate, the geometric con- 
tributions in the evolution will show dissipation-related 
aspects distinct from the dynamical evolution ones. A 
promising approach to address this issue is to foucs on the 
evolution of single trajectories. In this way one can re- 
sort to the theory of geometric phases developed for pure 
states. Individual trajectories can be realized condition- 
ally to the occurrence of a suitable chain of detections. 
Notice that this approach has been recently employed to 
calculate the geometric phase of open systems [lJi |2(| . 

A single trajectory in the quantum jump model is spec- 
ified, for a continuously measured system, as a chain of 
states obtained by the action of a specific sequence of 
detected jumps. The state evolution conditioned to the 
detection of no jumps, i.e., the no-jump trajectory, can 
be described by the Schrodinger equation 



i^m))=m)m)), 



(2) 



with the non-Hermitian effective Hamiltonian given by 



(3) 
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Here, the last term contained in Eq. accounts for 
the dissipative interaction. The state evolution indi- 
cated above therefore differs from the conventional uni- 
tary dissipation- free process. This deviation then results 
in an inevitable distortion for the desired gate operations. 

Let us first examine the Berry phase shift achieved 
through the double-loop refocusing scheme e.g., in an 
NMR system A spin half nucleus is assumed to 

be subjected to a cyclic field evolving adiabatically with 
cone angle 9 and its Hamiltonian is written as H§(t) = 
§ (cos 8cr z + sin9(cos^ta x + sin jta y )). The dynamical 
evolution will generate a dynamical phase <f>± — ±^QT 
as well as a geometric Berry phase = ±7r(l — cos 9). 
To remove the dynamical phase, the cyclic process is run 
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twice; the second loop is surrounded by a pair of hard 
7r pulses and it is driven in the opposite orientation of 
the first one. A pure geometric phase shift is then ob- 
tained. We consider now the system subjected to dissipa- 
tion with r = *Jkg-. The conditional no-jump evolution 
of the system is governed by the Schrodinger equation 
with an effective Hamiltonian ©. Instead of a real dy- 
namical phase originated from the unitary system, dis- 
sipation will give rise to a complex "dynamical phase" 
given by 

4> d ± =±~&T-^kT, ft = f>[sin 2 + (cos0-^) 2 ] 1/2 - 

2 2 ( 4 ) 

The imaginary part of 0± describes the decay of the tra- 
jectory probability a different decay rate is obtained for 
different basis states. This latter fact implies a distortion 
of the evolution respect to the dissipation- free process. 

More specifically for the desired geometric action, 
equation |0J reveals that the dissipation-induced distor- 
tion compared with <p 9 is in an order of k/7, which for 
an adiabatic process with small 7, can be rather severe. 
Nevertheless this distortion can be avoided in the whole 
procedure since, by removing the dynamical phase, the 
rcfocusing evolution with the opposite direction will re- 
move the inhomogeneity of dissipation as well. The over- 
all process gives merely origin to a global suppression of 
the trajectory probability with Yd = e~2 KT . 

The cancellation of the unwanted non-unitary effects 
achieved by the double loop technique is analog to a sort 
of noise decoupling by symmetrization via o^-pulses 0. 
Indeed the kinematical symmetry Hq '(t) = Hq(T — t) 
of the control process, along with the exchange of the 
computational basis states over the second loop leads to 
a homogeneous, i.e., global action of the dissipation that 
can then be neglected. 

It should be noted that the effect of dissipation on the 
dynamical phase discussed above is not the only influence 
that dissipation has on the quantum dynamics of the sys- 
tem. Indeed, the non-unitary character of the evolution 
will affect the geometrical phase too. One is then lead to 
consider a complex Berry phase describing the dissipative 
effects embedded inside the geometric action, 



±tt[i - {n/n)(cose- 
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In view of the sensitivity of the geometric phase under 
the reversal of the orientation of the second control loop 
the dissipative effects in Eq. (JSJ cannot be removed via 
the refocusing process. On the other hand, in case of 
low dissipation, equation (J5J reveals that the resulting 
error is of order ac/O and thus can be made negligible by 
making Q large enough. 

Dissipation in holonomic quantum computation. We 
turn now to consider another well-established approach, 
namely, the so-called holonomic quantum computation 
(HQC) 7]. In HQC information is encoded in a degen- 
erate eigenspace of the governing Hamiltonian. Pure ge- 
ometric operations, including the Abelian phase factor 



and non- Abelian transformation, are naturally achieved 
by adiabatic evolutions. Dynamical phase induced herein 
is clearly global due to the inherent structural symmetry 
(degeneracy) of the system levels. Dissipation occurring 
in this system can be divided into two categories: (i) 
dissipation preserving the degeneracy structure and (ii) 
inhomogeneous dissipation lifting the degeneracy. It can 
be easily seen that in case (i) of homogeneous dissipation, 
only global decay can be induced. Hence the distortion 
to the evolution can be avoided in the no-jump trajectory 
|23| . In what follows we shall focus on the general case 
of inhomogeneous dissipation, which will affect the cyclic 
evolution of the system in a more severe manner. 

Consider as an illustration the standard optical system 
|ToL IT^ I in which a qubit is encoded in a multi-level A- 
type trapped ion or a similar cavity atom. The ground 
(or metastable) levels \gt) (i — I, • • • ,n) are highly de- 
generated and each couples to the excited state |e) in a 
tunable manner. Let the states \gi) and 1 52) correspond 
to the computational basis |0) and |I), respectively. In 
order to obtain the gate operation e^ 9 ' 1 ^ 1 ', the system 
is described by the periodic Hamiltonian 

(6) 

where o e i := |e)(^| (i = 1,2,3), and the parame- 
ters 9, ip evolve adiabatically. In the dissipation-free 
case the dark state of the system given by \D(9,ip)) — 
cos 6*152)— sin ^e 11 * 3 1 (73), after aloop in the parameter space 
acquires a net Berry phase <f> 9 which depends only on the 
solid angle of the loop. Suppose now that the level \g^) is 
metastable and suffers damping with T^T — K3CT33. One 
has that {D{8, ip)\K,^azs\D{0, tp)) = K 3 sin 2 0. The dissi- 
pative dynamics lifts the degeneracy between the dark 
state and the ground state \gi). The evolution of the no- 
jump trajectory is described by the following non-unitary 
transformation 



W)) 
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where the operator uq is given by 



^+i0 9 )|i)<i 



(7) 



(8) 



with cf> d = —7} Jq K 3 sin 2 9dt. In detail, for the typical 
process in which ip — jt rotates uniformly, one has ex- 
plicitly tfi 9 — 47rsin 2 6 1 and <ft d = — tt(k 3 /j) sin 2 9. Dissi- 
pation induced distortion in this process scales as K3/7, 
which turns out to be unfavorable since for the adiabatic 
processing, the damping rate K3 may be comparable to 
the parameter frequency 7. Note also that in the above 
calculation, a low damping rate has been assumed: we 
have ignored the influence of dissipation on the geomet- 
ric phase 4> 9 and the deviation to the instantaneous dark 
state. Both of these effects would cause an error to the 
gate operation scaling as ^3/^. This latter ratio we as- 
sume can be made negligible. 
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Let us investigate further the influence of dissipation 
on the other single-qubit holonomic gate studied in [Io| . 
We consider the adiabatic evolution generated by the 
Hamiltonian 

H(t) = flsmO[cos(p(ai e + cr e i) + sillip((T2e + &e2)] 

+Qcos9(a 3e + a e3 ). (9) 

The system admits two dark states \Di) = 
cos0(cos<y9|gi) + sin(/?|g 2 )) — sin ^|_g 3 ) and \D 2 ) = 
cos ^152) — sint^lgi). A suitable cyclic evolution in the 
dissipation-free process will generate the non-Abelian 
holonomic gate operation e 1 ^ ay ; the parameter </> 9 
describes again the solid angle swept by the executed 
loop. We assume the same dissipation of the metastable 
level I (73). It should be noted that the dissipation 
generated dynamics for the two dark states does not 
commute with the non-Abelian connection. In view that 



(Di \n 3 a 33 \Dj) = K 3 sin for i 



J 



1 and null for 



other choices of i and j, a perturbation calculation gives 
directly that 



uc 



a d |O)<O|+# s 



(10) 



where (f> d — k 3 sin 2 9dt. For the specified pro- 

cess with ip — yt, one has <fi 9 — 2-7rcos6' and cj> d — 
— ""(^3/7) sin 2 9. Similar to the former case, the distor- 
tion induced by inhomogeneous dissipation to this gate 
is also of order K3/7. 

Generalized refocusing scheme for holonomic quantum 
computation. We have shown in the above the influ- 
ence of dissipation on two different schemes for geomet- 
ric quantum computation. At this stage, an interesting 
question arises: Is the refocusing strategy described in 
the former scheme extendable to the case of non-Abelian 
holonomic operations? In what follows we will explic- 
itly address that question for the kind of optical system 
proposed in |ldj| . The answer will turn out to be "yes" . 

It is worthwhile to stress the role of symmetry in the 
dynamical evolutions of above systems. The kinematical 
symmetry of refocusing scheme, i.e., H c (t) = H £ (T — t), 
helps in quenching the dissipation-induced distortion. 
This is due to the fact that refocusing removes the dy- 
namical contributions from the evolution and these latter 
are the ones which are more severely affected by dissipa- 
tion. On the other hand, the symmetry associated with 
the level degeneracy in the holonomic scheme is lifted 
by the dissipative interaction. This in turn gives rise to 
an inhomogeneous decay of basis states and to a corre- 
sponding distortion of the holonomic gates. The scheme 
we are going to discuss in the following basically combines 
the kinematical symmetry of the double-loop refocusing 
scheme with a technique for restoring the dark-state sym- 
metry and therefore homogeneity of the decay process. 

In brief the novel refocusing scheme for holonomic 
quantum computation that we are going to analyze, can 
be described as a double-loop configuration in which 
the additional loop refocuses the first one in an oppo- 



site direction. As one may see, the opposite loop exe- 
cution would induce the inverse operation hence cancel 
the desired transformation, but this can be side stepped 
through interchanging the roles of the two computational 
bases in the twice loop evolutions. 

Let us start by considering in a detailed fashion the 
single-qubit gate e 1 ^ 9 ' 1 ^ 1 '; in the new scheme this gate 
can be achieved simply by performing the loop evolution 
twice, the first one as in (j^Jl followed by a second one 
described by 

H c (9,p) = flsm9(a le + a el ) + flcos9(a 3e e iv + a e3 e~ iv ), 

(11) 

where the parameters 9 and p refocus the first loop but 
in the opposite direction. Note in the evolution of Eq. 

the level \gi) was employed to replace |<?2) in the 
first loop of Eq. . We have interchanged the roles of 
the two basis states \gi) and \g2) in order to make the 
overall dynamical evolution symmetric with respect the 
exchange of states of the computational basis. This sort 
of symmetrization plays the role of the two hard tt pulses 
used in the former NMR scheme to flip the basis states. 

The dissipative dynamics on the second loop l|llf> will 
generate the transformation with 



„(^-i^)|0><0| 



(12) 



Thus the total two-loop process will be described by the 
following overall transformation 



Udl 



e d e -^ 9 e i2^|l)(l| 



(13) 



where e~ 1 ^ 3 denotes an irrelevant global phase. From 
this equation it can be clearly seen, in view of the sym- 
metric role played by |0) and |1), that dissipation gives 

rise just to a global factor affecting just the no-jump 
trajectory probability; the former distortion in equation 
(JSJ) has been suppressed. 

We now discuss the extension of the refocusing strategy 
to the second single-qubit gate considered in 01 , i- e -> 
e j0 9 <T H . fois extension turns out to be more difficult than 
the previous one. Indeed, one may naively expect that 
the refocusing loop is described by 

H<j{9,p) — r2sin0[cosy>(<T2e + °~ei) + sinp(ai e + cr e i)] 
+Q cos 9{a 3e + cr e3 ) (14) 

with the parameter 9 and <p evolving in the opposite di- 
rection of the first loop 0. However, it is immediate 
to verify that the evolution associated to the loop l|14fl . 
namely the transformation uq = e^ d l 1 ^ 1 l +l ^ 9cr « , does not 
commute with that of Eq. (|10fl . In other terms the non- 
commutativity of the two loop transformations destroys 
the desired symmetry aimed by the refocusing process. 
As a result such a double-loop evolution does not of- 
fer any clear advantage as far as resilience against dis- 
sipation is concerned. One can overcome this obstacle 
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by performing the two opposite controlled evolutions in 
sense simultaneously. As the superposed process shall 
lead overlapped interactions to the level I53), an extra 
level will be necessary to pin it down. Explicitly, the 
whole process of the dual-superposed loop evolution can 
be described by the following periodic Hamiltonian 

Hdi = r2sin0[(cost^ — sint/?)cri e + (sinyj + cos(p)(T2 e ] 
+Q cos 9(<j3 e + cr 4e ) + h.c. (15) 

Here the level \gs) and the extra level \g±) are supposed 
to suffer the homogeneous damping T'T = K3(o33 + 0-44). 
As the parameters 9 and ip evolve adiabatically, the two 
dark states of the newly established Hamiltonian 1)15(1. 

\Df) = cos0(cosy%i) +smip\g-2)) -sin% 3 ), 
\Df) = cos 8 (cos ip\g 2 ) - sux(p\gi}) - sin% 4 ), (16) 

shall generate a non-Abelian holonomy and it can imple- 
ment the identical gate operation e l ^ s ' Jy with <j) 9 denot- 
ing the loop-related solid angle. The advantage of the 
present system is obvious: the yielded new dark states 
possess completely symmetric structure which is thus ca- 
pable to remove the inhomogeneous effect of dissipation 
on the computational space. Simply the dissipation can 
merely lead to a global decay and for the no-jump tra- 
jectory the transformation is described by 




It should be pointed out that the loop superposition 
strategy is also realizable for the previously discussed 



single-qubit holonomic gate, which thus provides an al- 
ternative way to implement the gate e 1 ^ WW, We 
stress here that the homogeneity achieved in the loop- 
superposed processes is restricted on the two-dimensional 
subspace spanned by the computational basis. 

In this paper we analyzed the robustness of geometric 
quantum gates, both Abelian and non-Abelian, against 
quantum dissipation. We have adopted the quantum 
jump approach to describe dissipative dynamics and we 
focused on the no-jump trajectory. It has been shown 
that geometric gates realized by means of a refocused 
double-loop scheme possess a certain resilience against 
dissipation-induced distortion. This robustness is due to 
the underlying symmetry of the scheme with respect to 
the exchange of computational basis states. The cru- 
cial point is that the non-unitarity affects mostly the dy- 
namical phases; these latter can be removed by a suit- 
able control process which leaves the geometrical part 
of the evolution untouched. The power of this refocus- 
ing strategy has been demonstrated by extending it to 
non-Abelian holonomic gates. For the non-Abelian case, 
we have shown how the same sort of symmetrization can 
be achieved via a loop superposition strategy. We con- 
clude by noticing the central role of symmetry in a variety 
of schemes of quantum information processing (see e.g., 
Ref. [24| and references therein). The present work re- 
veals one more intriguing way in which symmetry might 
help in our struggle against decoherence. 
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[1] P.W. Shor, Phys. Rev. A 52, 2493 (1995); A.M Steane, 
Phys. Rev. Lett. 77, 793 (1997). 

[2] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 

[3] D. Gottesman, Phys. Rev. A 54, 1862 (1996). 

[4] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 
(1997); Mod. Phys. Lett. B 11, 1085 (1997). 

[5] D.A. Lidar, I.L. Chuang, and K.B. Whaley, Phys. Rev. 
Lett. 81, 2594 (1998). 

[6] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998); L. 
Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 2417 
(1999); P. Zanardi, Phys. Lett. A 258, 77 (1999). 

[7] P. Zanardi and M. Rasetti, Phys. Lett. A 264, 94 (1999). 

[8] J. A. Jones, V. Vedral, A. Ekert and G. Castagnoli, Na- 
ture (London) 403, 869 (2000). 

[9] G. Falci, R. Fazio, G.M. Palma, J. Siewert and V. Vedral, 
Nature (London) 407, 355 (2000). 
[10] L.M. Duan, J.I. Cirac and P. Zoller, Science 292, 1695 
(2001). 

[11] J. Pachos, P. Zanardi, and M. Rasetti, Phys. Rev. A 61, 

010305(R) (2002). 
[12] A. Recati, T. Calarco, P. Zanardi, J.I. Cirac and P. Zoller, 

Phys. Rev. A 66, 0302309 (2002). 
[13] L.X. Cen, X.Q. Li, Y.J. Yan, H.Z. Zheng, and S.J. Wang, 



Phys. Rev. Lett. 90, 147902 (2003). 
[14] M.V. Berry, Proc. R. Soc. London, Ser. A 392, 45 (1984); 

B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 
[15] F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984). 
[16] A. Nazir, TP. Spiller, and W.J. Munro, Phys. Rev. A 

65, 042303 (2002). 
[17] A. Blais and A.-M. S. Tremblay, Phys. Rev. A 67, 012308 

(2003). 

[18] G. DeChiara and G.M. Palma, Phys. Rev. Lett. 91, 
090404 (2003). 

[19] A. Carollo, I. Fuentes-Guridi, M.F. Santos, and V. Ve- 
dral, Phys. Rev. Lett. 90, 160402 (2003). 

[20] A. Carollo, I. Fuentes-Guridi, M.F. Santos, and V. Ve- 
dral, Phys. Rev. Lett. 92, 020402 (2004). 

[21] I. Fuentes-Guridi, F. Girelli, and E. Livine, 
|quant-ph/0311164| 

[22] M.B. Plenio and P.L. Knight, Rev. Mod. Phys. 70. 101 
(1998). 

[23] In case of homogeneous dissipation, the trajectories with 
jumps occurring have also been discussed in Ref. [2l|. 

[24] P. Zanardi, and S. Lloyd, Phys. Rev. Lett. 90, 067902 
(2003). 



